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■ Abstract The Degasperis - Procesi (DP) equation describing the propagation of 

shallow water waves contains a physical parameter u, and it is well-known that the 
DP equation admits solitary waves with a peaked crest when u = 0. In this article, we 
illustrate, for the first time, that the DP equation admits peaked solitary waves even 
when oj ^ 0. This is helpful to enrich our knowledge and deepen our understandings 
about peaked solitary waves of the DP equation. 
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in '. 1 Introduction 
in 

£N) ■ The solitary surface wave was first discovered by Russell [17] in 1834. Since then, 
many models for shallow water waves have been developed, such as the Boussinesq 
equation [3J, the Korteweg & de Vries (KdV) equation [5] and the Benjamin-Bona- 
^ Mahony (BBM) equation [TJ, and so on. All of them admit solitary waves with 

smooth crest. In contrast to the KdV, Boussinesq and BBM equations, the celebrated 
Camassa-Holm (CH) equation [4] 

Ut + 2uu x — u xx t + 3uu x = 2u x u xx + uu xxx (1) 

can model both phenomena of soliton interaction and wave breaking (see [5]), where 
u(x,t) denotes the wave elevation, x,t are the temporal and spatial variables, u is a 
constant related to the critical shallow water wave speed, and the subscript denotes 
the partial differentiation, respectively. The CH equation ([T]) is integrable and bi- 
Hamiltonian. Especially, when u — 0, the CH equation (TjQ) has the peaked solitary 
wave [1] 

u(x, t) = cexp(— \x — ct\). 
Degasperis and Procesi [7J investigated a generalized equation 

u t + 2uu x - u xxt + {(3 + l)uu x = (3u x u xx + uu xxx (2) 



2 



and found that it is completely integrable only when (3 = 2, corresponding to the CH 
equation JT]), and (3 = 3, corresponding to the so-called Degasperis - Procesi (DP) 
equation 

Ut + 2uu x - u xxt + Auu x = 3u x ii xx + uu xxx , (3) 

respectively. The DP equation ([3]) can model the propagation of shallow water waves 
with small amplitude and large wave-length. It is interesting that, like the CH equa- 
tion (Q]), the DP equation (j3J) also admits solitary waves with a peaked crest when 
u = 0. 

Currently, Liao [13J gained the closed-form solutions of the peaked solitary waves 
of the KdV equation, the modified KdV equation, the Boussinesq equation and the 
BBM equation. Besides, by means of the Mathematica package BVPh 1.0 that is 
based on the homotopy analysis method (HAM) [5HT2]. a powerful analytic method 
for highly nonlinear equations, Liao [H] found that the CH equation ([1]) also admits 
peaked solitary waves even in the case of w ^ 0. All of these indicate that most of 
mainstream models of shallow water waves admit peaked solitary waves. 

In this article, using the HAM-based Mathematica package BVPh 1.0 for nonlinear 
boundary value problems in a similar way, we illustrate that the DP equation (|3]) also 
admits peaked solitary waves even when u ^ 0. In addition, we further illustrate that 
Eq. 02]) for arbitrary constant (3 and the more generalized equation 

u t + 2uu x - u xxt + auu x = (3u x u xx + uu xxx (4) 

for arbitrary constants a and (3 also admit peaked solitary waves even when u ^ 0, 
although they are not integrable in general cases. 



2 Peaked solitary waves of Eq. (EE) in the case of 

u ^ 

Let us first consider the propagating solitary waves of the DP equation ([3]) with 
permanent form, corresponding to a = 4 and = 3 of Eq. (T4]). To avoid the repeat, 
Eq. 01]) is used to describe the basic ideas of our analytic approach. Writing £ = x — ct 
and u>(£) = c u(x,t), the original equation (TjJ becomes 

w '» _ ( i _ _^ ) w > + aww > = p w > w » + ww >"^ (5) 



c / 

subject to the boundary conditions 

w ->• 0, w' -> 0, w" -> 0, as |£| -> +oo, (6) 
where the prime denotes the differentiation. 
The corresponding linearized equation 

w'" -(l-^pjw' = (7) 
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has the peaked solitary wave 



w(0 = A exp(-^|), ix = y 1 - ^, (8) 

provided u < c/2, where A = w(0). The solution of Eq. (jSJ) can be expressed in the 
form 



'(0 = *>o(0 + 5>m(0- ( 9 ) 



n=l 

where a n is a consant to be determined. Using the homotopy analysis method (HAM) 
[91412]. an analytic technique for highly nonlinear differential equations, we gain the 
series solution 

+oo 
m=0 

Here 

w (0 = Aexp(-fi\£\) 
is the initial guess, w m for m > 1 is governed by 

£ Kn(£) - Xm W m _i(0] = C 5m~l(0, ( 10 ) 

subject to the boundary condition 

w m (0) = 0, w m ->■ 0, as |f| -> +oo, (11) 

where 

is an auxiliary linear operator, and 

n 

i=o 

J 1 when n > 1, „s 
Xn ~ \ otherwise. ^ J 

Note that Co 7^ is an auxiliary parameter, called the convergence-control parameter, 
which provides us a convenient way to guarantee the convergence of the approximation 
series. For details, please refer to Liao [U[TT1[T2]. In fact, directly using the HAM- 
based Mathematica package BVPh 1.0 (see Part II of Liao's book [12]) for nonlinear 
boundary- value/eigenvalue problems, it is straightforward to gain high-order analytic 
approximations of (jSJ) and (JfJJ). For details, please refer to the Appendix. 

The accuracy of the mth-order approximation of w (£) is defined by the averaged 
residual square of the governing equation (0) in the domain £ e [0, a]: 



£ _ 

a 



"IE 

3=0 



d£, (14) 



4 



Table 1: u'(0+) of the mth-order analytic approximations and the residual squares 
of the DP equation (j3J) in the case of c = 1,cj — 1/4, and A = ±1/10 given by the 
HAM-based Mathematica package BVPh 1.0 with the convergence-control parameter 
c = -1. 



m 


A = 1/10 


A = -1/10 


Ti'(Q+) 


P 


u'(0+) 


P 


2 


-0.065250 


3.8 E-9 


0.074675 


4.8E-9 


4 


-0.065110 


2.0 E-12 


0.074779 


1.6E-12 


6 


-0.065105 


9.2 E-15 


0.074782 


5.9E-15 


8 


-0.065105 


2.4 E-17 


0.074782 


1.6E-17 


10 


-0.065105 


5.2 E-20 


0.074782 


3.5E-20 



where 

Afw = w'" — fi 2 w' + a ww' — (3 w'w" — ww'". 

We choose a = 10 in this article, because the wave elevation decays exponentially. 

For simplicity, we investigate the case c = 1 in this article. Without loss of 
generality, we first consider the case of u = 1/4 and A = ±1/10 for the DP equation 
([3]). As shown in Tabled], the averaged residual squares of the lOth-order analytic 
approximations decrease to 5.2xl0~ 20 in the case of A = 1/10 and to 3.5xl0 -20 in 
the case of A = —1/10, respectively. Besides, the corresponding values of u'(0 + ) (the 
limit is taken as £ — > from the right) quickly converge to -0.065105 and 0.074782, 
respectively. All of these clearly indicate that the series given by the Mathematica 
package BVPh 1.0 converge to the solution of the DP equation ([3]). The solutions 
for larger values of \A\ can be gained in a similar way. All of these solitary waves 
have a peaked crest, as shown in Fig. [TJ These clearly illustrate that, like the CH 
equation ([TJ), the DP equation (j3J) also admits peaked solitary waves even when u^0. 
Besides, like the CH equation ([TJ, the phase speed of the peaked solitary waves of 
the DP equation (J3J) has nothing to do with the wave amplitude A. This is quite 
interesting. 

As pointed out by Degasperis and Procesi [TJ, the generalized equation ([2J) is 
completely integrable only when p = 2 or p = 3, corresponding to the CH equation 
([TJ) or the DP equation (J3J), respectively. Does Eq. <$Z§ admit peaked solitary waves 
when p ^ 2 and p ^ 3 ? 

Without loss of generality, let us consider the case ofc = l,a; = l/4, J 5 = 5/2 and 
A = ±1/10 for Eq. (T5J), corresponding to a = 7/2,(3 = 5/2 for Eq. (T4J). Similarly, by 
means of the HAM-based Mathematica package BVPh 1.0, we gain the convergent 
approximations with cq = — 1. As shown in Table [2j the averaged residual squares 
of the lOth-order analytic approximations decrease to 4.4xl0 -21 in the case of A = 
1/10 and to 3.2xl0 -21 in the case of A = —1/10, respectively. Similarly, using the 
Mathematica package BVPh 1.0, we gain the convergent approximations for larger 
values of A, as shown in Fig. [21 Note that all of these solitary waves have a peaked 
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Figure 1: The peaked solitary waves u(x,t) of the DP equation (J3J) in the case of 
u = 1/4 and c = 1. Solid line: A = ±3/10; Dashed line: A = ±1/5; Dash-dotted 
line: A = ±1/10. 

crest! Thus, the generalized equation (J2J) in the case of w / also admits peaked 
solitary waves, even when (3 ^ 2 and /3 ^ 3 so that it is not integrable! 

How about the more generalized equation 

Without loss of generality, let us consider the generalized equation (j3J) in the case 
of c = l,u = 1/4, a = 11/3, /5 = 5/2 and A = ±1/10. Similarly, by means of the 
HAM-based Mathematica package BVPh 1.0, we gain the convergent approximations 
with Co = — 1. As shown in Table El the averaged residual squares of the 10th- 
order analytic approximations decrease to 1.1 xlO -20 in the case of A = 1/10 and to 
7.8xl0 -21 in the case of A = —1/10, respectively. Similarly, using the Mathematica 
package BVPh 1.0, we gain the convergent approximations for larger values of \A\, as 
shown in Fig. [31 It should be emphasized that all of these solitary waves have a peaked 
crest. These illustrate that the generalized equation (J4J) also admits peaked solitary 
waves for arbitrary constants a and (3 even when w^0, although the corresponding 
equation is not integrable in general cases. Therefore, the existence of peaked solitary 
waves might be a common property of the generalized equation (T4J). 



As pointed out by Constantin and Molinet [6], all of these peaked solitary waves 
should be understood mathematically as weak solutions of the generalized equation 
(jlj). However, physically, this kind of discontinuity of wave elevation widely appears 
in fluid mechanics, such as dam break [19] in hydrodynamics and shock waves in 
aerodynamics, which have clear physical meanings. Note also that such kind of dis- 
continuous problems belong to the so-called Riemann problem [2|[T6 | [T8 | [T9] , a classical 
field of fluid mechanics. 
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Table 2: u'(0+) of the mth-order analytic approximations and the residual squares 
of Eq. P) in the case of c = l,u = 1/4, A = ±1/10 and = 5/2 given by the 
HAM-based Mathematica package BVPh 1.0 with the convergence-control parameter 
c = -1. 



m 


A = 1/10 


A = -1/10 


u'(0+) 


£ 


u'(0+) 


£ 


2 


-0.066002 


2.1 E-9 


0.074251 


2.6E-9 


4 


-0.065904 


5.7 E-13 


0.074323 


5.3E-13 


6 


-0.065901 


1.6 E-15 


0.074325 


1.1E-15 


8 


-0.065901 


3.0 E-18 


0.074325 


2.2E-18 


10 


-0.065901 


4.4 E-21 


0.074325 


3.2E-21 



Table 3: u'(0 + ) of the mth-order analytic approximations and the residual squares 
of the generalized equation (@J in the case of c = l,u — 1/4, A = ±1/10, a = 11/3 
and (3 = 5/2 given by the HAM-based Mathematica package BVPh 1.0 with the 
convergence-control parameter cq = —1. 



m 


A = 1/10 


A = -1/10 


u'(0+) 


£ 


u'(0+) 


£ 


2 


-0.065541 


2.3 E-9 


0.074576 


2.9E-9 


4 


-0.065428 


9.8 E-13 


0.074658 


7.6E-13 


6 


-0.065425 


3.3 E-15 


0.074660 


2.2E-15 


8 


-0.065425 


6.6 E-18 


0.074660 


4.6E-18 


10 


-0.065425 


I A E-20 


0.074660 


7.8E-21 



3 Concluding remarks 

In this paper, using the HAM-based Mathematica package BVPh 1.0 as a useful 
tool, we illustrate that, like the CH equation ([TJ, the DP equation ([3]) also admits 
peaked solitary waves even when u ^ 0. Furthermore, we illustrate in a similar way 
that the generalized equations (EJ and PJ also admit peaked solitary waves even 
when u ^ 0, and the corresponding equation is not integrable! As pointed out by 
Liao [13] , nearly all mainstream models of shallow water waves admit peaked solitary 
waves. This work further confirms this view-point. This is mainly because the fully 
nonlinear water wave equations might admit peaked solitary waves, as pointed out 
by Liao [15] . Therefore, mathematically speaking, peaked solitary waves widely exist 
and might be a common property of models for water waves, although it should be 
carefully investigated whether or not such kind of peaked solitary waves indeed could 
be regarded as weak solutions and besides have physical meanings. 



Figure 2: The peaked solitary waves u(x, t) of Eq. <^ in the case of uj — 1/4, c = 1 
and = 5/2, corresponding to a = 7/2,(3 = 5/2 of Eq. (0}. Solid line: A = ±3/10; 
Dashed line: A = ±1/5; Dash-dotted line: A = ±1/10. 

Acknowledgement 



This work is partly supported by the State Key Lab of Ocean Engineering (Approval 
No. GKZD010056-6) and the National Natural Science Foundation of China. 



8 




Figure 3: The peaked solitary waves u(x, t) of (jlj) in the case of u = 1/4, c = 1, 
a = 11/3 and p = 5/2. Solid line: A = ±3/10; Dashed line: A = ±1/5; Dash-dotted 
line: A = ±1/10. 
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Appendix 

The use of HAM-based Mathematica package BVPh 1.0 



The Mathematica package BVPh 1.0 for nonlinear boundary- value/eigenvalue 
problems is developed and issued by Liao [12] (Part II), which is based on the HAM 
and is free available online. Using the BVPh 1.0, it is straightforward to gain the 
analytic approximations of the peaked solitary waves of Eqs. flSJ) and (E]) for given 
w(0) = A. Here, we briefly describe how to do it in case ofc = l,w = l/4, a = 4,/3 = 3 
and A = 1/10, corresponding to the DP equation (j3]). 



1. First, download the BVPh 1.0 (the code file is named by BVPh_L0.txt) online 



( http://numericaltank.sjtu.edu.cn/BVPh.htm ) and save it in a directory such 



as C:/math/DP as an example. 

2. Then, run the computer algebra system Mathematica, and type the following 
command one by one: 



SetDirectory ["C:\math\DP"] ; 



«InputDP.txt 



The file named InputDP.txt contains the following Mathematica commands and 
necessary definitions for BVPh 1.0: 



(* Install the BVPh 1.0 *) 
«BVPhl_0.txt; 



(* Define the physical and control parameters *) 
TypeEQ = 1 ; 

ApproxQ = ; 

ErrReq = 10" (-30); 

zRintegral = 10; 



(* Define the governing equation *) 

ALPHA = 4 ; 

BETA = 3; 

mu2 = l-2*omega/c ; 

f [z_,u_,lambda_] := D [u, {z ,3}] -mu2*D [u,z] \ 

+ ALPHA*u*D[u,z] - BETA*D[u,z]*D[u,{z,2>] - u*D[u,{z,3}] ; 
If [BETA == 2 && ALPHA == 3, Print ["CH equation"]]; 
If [BETA == 3 && ALPHA == 4, Print ["DP equation"]]; 



(* Define Boundary conditions *) 
zR = Infinity; 
DrderEQ = 3; 

BC[l,z_,u_,lambda_] := Limit[u-A, z -> ] ; 
BC[2,z_,u_,lambda_] := Limit [u, z -> zR ] ; 
BC[3,z_,u_,lambda_] := Limit [D [u,z] , z -> zR ] ; 

(* Define initial guess *) 

mu = Sqrt [mu2] ; 

u [0] = A*Exp [-mu*z] ; 

(* Define output term *) 

output [z_,u_,k_] := Print ["output = " , D [u [k] , z] /. z->0//N] ; 

(* Defines the auxiliary linear operator *) 
L[u_] := D[u,{z,3}] - mu2 * D[u,z]; 

(* Print input and control parameters *) 
Print Input [u [z] ] ; 

(* Set convergence-control parameter cO and physical parameters 
cO = -1 ; 
A = 1/10; 
omega = 1/4; 
c = 1; 

Print [" cO = ",c0, " omega = ", omega, " c = ", c, " A 



(* Gain up to lOth-order HAM approxiamtion *) 
BVPh[l,10] ; 

(* Get results in the whole domain *) 
For[k=0,k<=10,k++,W[k] = U[k] /. z-> Abs[x]]; 

(* Show the 5th and lOth-order approximation *) 

Plot [{W [5] , W [10] } , {x , - 10 , 10} , PlotRange->{Min [A , 0] , Max [A , 0] }] 



